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We carry out a systematic study of the dispersion relation for linear electrostatic waves in an 
arbitrarily degenerate quantum electron plasma. We solve for the complex frequency spectrum for 
arbitrary values of wavenumber k and level of degeneracy /r. Our finding is that for large k and 
high fj, the real part of the frequency oJr grows linearly with k and scales with only because of 
the scaling of the Fermi energy. In this regime the relative Landau damping rate 'y/uur becomes 
independent of k and varies inversly with fj,. Thus, damping is weak but hnite at moderate levels of 
degeneracy for short wavelengths. 


Introduction. Although the field of plasma physics 
has traditionally been concerned with systems of clas¬ 
sical particles, the treatment of collective interactions 
of charged quantum particles has been of interest for 
many years [iH^. More recently, interest in quantum 
plasma physics has increased due to the possible im¬ 
portance of collective plasma effects in microelectronic 
systems (a; the experimental realization of warm, dense 
matter [ 6 |; and in the description of astrophysical com¬ 
pact objects Q. 

Generally, quantum effects are expected to be sig¬ 
nificant for physics at the smallest scales. For plasma 
physics, this means systems in which microscopic scales 
such as the Debye length \d = {T/Airne^) and gyro- 
radius Pf, = meCv±/eB are on the order of or smaller 
than the de Broglie wavelength of the particles, X^b = 
h/ {fl-nmeT) ' . The proper treatment of such micro¬ 
scopic scales requires kinetic theory. For this reason, de¬ 
veloping a kinetic theory of quantum plasma physics is 
of special interest. 

One classic textbook problem of classical kinetic 
plasma physics is the propagation, dispersion, and col¬ 
lisionless Landau damping of linear electrostatic elec¬ 
tron waves in unmagnetized plasmas [ 8 |. This problem 
has proven useful in illustrating the effects of quantum 
physics on plasma phenomena as well [i[i3- Landau 
damping in the quantum case is also of interest as it 
is potentially observable in experiments involving nano¬ 
scale electronics, plasmonic devices and warm dense mat¬ 
ter [^, ITll-fl3l | and of importance in dense astrophysical 
plasmas |l4| . Plasmon dispersion and damping has been 
observed in solid-state plasmas [iM3- On the theo¬ 
retical front, the linear quantum longitudunal dispersion 
relation has been studied extensively in the literature in 
a variety of limiting cases, including: the c omp letely de- 
gen erate limit with phase velocities above |lOl | and near 
(l§ . [l^ the Fermi velocity; at arbitrary phase velocities 
in one-dimensional plasmas for several temperatures ; 
and the high phase velocity limit for arbitrarily degener- 
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ate quantum plasmas including the ion-acoustic 

mode [13 . In addition, there have recently been stud¬ 
ies of the onset of nonlinearity and particle trapping in 
quantum Landau damping |23l |24| . However, a system¬ 
atic examination of the linear dispersion relation for all 
wavelengths in arbitrarily degenerate plasmas has so far 
not appeared in the literature. There is, for example, no 
study of short-wavelength wave propagation or damp¬ 
ing in moderately degenerate plasmas. For this reason, 
a thorough study of the dispersion of linear electrostatic 
waves of any wavelength in both partially and completely 
degenerate quantum plasmas is called for. 

In this paper, we study linear one-dimensional (ID) 
longitudinal waves in an arbitrarily degenerate, unmag¬ 
netized electron gas, with a neutralizing uniform ion 
background. We map out numerically the complex dis¬ 
persion relation over a broad range of wavenumbers k 
and levels of degeneracy p/T, p being the chemical po¬ 
tential and T the temperature. We then use this map 
of the complex frequency w = — 17 to analyze the 

key differences between non-degenerate and degenerate 
cases and discuss the possibility of Landau damping in 
the degenerate limit. We find that for long wavelengths 
LOr is approximately constant and the damping rate 7 is 
exponentially suppressed for all values of degeneracy. At 
short wavelengths however, both ojr and 7 grow linearly 
with k so that the relative Landau damping rate 'y/ujr 
becomes independent of k. The short-wavelength rel¬ 
ative damping rate exceeds unity in the classical case, 
but decreases as {p/T)~^ as p —>■ 00 . This implies that 
it may be important to account for both the presence 
of electrostatic waves and their damping in the descrip¬ 
tion of a real partially degenerate plasma, depending on 
the specific level of degeneracy. This additionally implies 
that electrostatic waves of arbitrarily short wavelength 
are able to propagate in a degenerate plasma. 

Theory of Linear Electrostatic Waves in a Quan¬ 
tum Plasma. Quantum kinetic plasma physics can be 
cast in a particularly convenient form by utilizing the 
phase-space formulation of quantum mechanics (^ . 
One defines the 1-body Wigner function F (x, p; f) = 
Wip) = J (x -k y |/3| X - y) [H, where p = 

I'*/’*) ('*/’/ density operator. The operation W, 

called the Wigner transformation, defines the quantum 
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phase space and F acts as a quasi-distribution function. 
For a system with a scalar potential V and Hamil- 
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tonian H = 127X1 + V, the Wigner function evolves 

according to the Moyal equation 
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where the subscripts indicate the function upon which 
a given derivative operates and the sine of differential 
operators is defined by its power series. This equa¬ 
tion becomes the Vlasov equation in the classical limit 
h ^ 0, with corrections of higher order in h acting as 
scattering/non-phase-space-volume-preserving terms. 

Using Eq. ([1]) linearized about a background equilib¬ 
rium distribution function Fq, one can obtain the longi¬ 
tudinal dielectric function [l|, 0, EUl, 
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where A = fik/rrie and ujp = (^Airne^/me) is the clas¬ 
sical electron plasma frequency. We integrate Eq. [5] 
over the directions perpendicular to the wave propaga¬ 
tion k, obtaining a ID problem with a reduced equilib¬ 
rium distribution function of the velocity parallel to k, 
/o (t'li) = / (Pv±Fo (v). The remaining integral is along 
a contour C following the real uy axis and diverting into 
the negative Im (w||) half-plane to encompass the singu¬ 
larity at U|| = uj/k. This diversion must be done in such 
a way that the contour also avoids any singularities in 
the function /q. 

The dielectric function in Eq. m differs from the 
classical case in two ways. First, the velocity deriva- 
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tive of the distribution function is replaced with a fi¬ 
nite difference in the plasmon momentum hk. Sec¬ 
ond, quantum mechanics can influence the background 
distribution Fq, which for arbitrarily degenerate plas¬ 
mas in thermal equilibrium is a Fermi-Dirac distribution 
Fo(v,/r) oc [ 1 -I-exp (v^/2meT — p/T)] The corre¬ 
sponding reduced ID distribution, defined by integrating 
over the perpendicular directions, is 

A («,, ^) = Li, (1 + _ , 3 ) 


where N (p) = [^ 7103/2 (-e'^/^)] ^ is the degeneracy- 
dependent normalization, thus normalizing /o to 1, and 
VT = (T/me)^^^ is the classical electron thermal veloc¬ 
ity. Here Li 3/2 is the poly logarithm function of order 3/2. 
The ratio p/T determines the level of degeneracy of the 
system: a Maxwellian is recovered in the limit p/T —>■ 

—00 and a distribution with Jq (y > vp) = 0 in the de- 

1 /2 

generate limit p/T —>■ -foo, where vp = {/2Ep/me) ' 
is the Fermi velocity with Ep = {/h?/2me) (2?™)^^^ the 
Fermi energy. 

By making a transformation to nondimensional vari¬ 
ables X = v\\/vt, k —>■ kXp, and w —>■ w/wp in Eq. ([2) 
along with Eq. we obtain 
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Im (w) < 0 


Here P is the principal value, ^ = exp(/x/T), and F[ = 
/‘^Xp = 87r-\/2Swp/T is the dimensionless quan¬ 
tum recoil parameter. 

In Eq. (|4]), there are three parameters: k, ^ (or /x), 
and F[. The important physical length and time scales 
are Xp and ■ Additionally, there exists a natural 
scale for wave phenomena given by the screening length 
A, = {v)/ujp where (v) = f |v| Fq (v,/x/T) is the 
average speed of particles in the plasma. The level of 


degeneracy determines (v) and thus the screening length 
is dependent upon degeneracy, 

A*=Li5/2(-e'^/^)Ai5, (5) 

which becomes the Debye length and Thomas-Fermi 
length Xtf = vp/ujp in the non-degenerate and degener¬ 
ate limits, respectively. The dependence of this screening 
length on the chemical potential is shown in the inset in 
Fig. [TJ All of the above parameters and scales depend 
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on n and T. 

The physical parameters of the system we describe 
are the density and temperature n and T, but these 
parameters are manifest in the dielectric function Eq. 
(4) through the dimensionless parameters ^ (n, T) /T and 
H {n, T). For simplicity we have limited this study to the 
case of small H = 10“^, varying only /r while keeping H 
fixed. This means that n and T cannot take arbitrary 
values. In future work it will be beneficial to map the 
dispersion relation for arbitrary values of n and T (and 
hence of p. and iJ), but for now we focus on the mathe¬ 
matical properties of the dispersion relation without con¬ 
sidering all possible physical parameter regimes. 

The quantum kinetic theory utilized in this paper 
is applicable for weakly-coupled non-relativistic plasmas 
without spin. The first condition requires a large num¬ 
ber of particles to be present within a screening radius 
n\^ ^ 1 , and the second condition requires that the 
thermal and Fermi velocities not approach the speed of 
light. In the zero temperature (totally degenerate) case 
this has the effect of setting minimum and maximum val¬ 
ues for the density, namely nX^p = y^ir^nap ^ 1 (where 
qb = lrne,e^ is the Bohr radius) and vp/c <C 1. Nu¬ 
merically this means 2.2 x 10^^cm“^ <C n <C 7 x 10^®cm^ 
The extension to relativistic and spin plasmas is an on¬ 
going effort. 

Results. We evaluate the dielectric function e from 
Eq. (|^ by direct numerical integration. Taking the usual 
logarithm branch cut in the negative horizontal direction, 
all integrals along the real axis are well-defined. We first 
solve e {oj, k, fi, H) = 0 for uj for single values of k, ^,/T 
and H in the classical and uj/k ^ vt limits using New¬ 
ton’s Method with a starting point given by an analytical 
approximation. We then solve for the entire classical dis¬ 
persion relation by iterating to larger values of k. We 
then iterate to larger values of /r/T using the results of 
the next-smallest degeneracy. In this way we map out oj 
in the entire (fc, /r) plane with fixed value oi H = 10“^. 

A complication arises due to the dependence of the 
screening length A* on /i, varying between Ad in the clas¬ 
sical case and Xtf (m) in the degenerate regime. We thus 
consider the dependence of ujr and 7 on /r and fcA, instead 
of kXo- This dependence is shown for several values of 
degeneracy in Fig. [T] In all cases, for small k, uJr is 
approximately cOp and weakly dependent on fc, and for 
fcA* > 1 grows linearly with k. The damping rate 7 is 
exponentially suppressed for small k at all degeneracies, 
and also grows linearly at large k. The effect of increas¬ 
ing degeneracy on ui (k) is to change the sharpness of the 
transition between flat and linear growth of ujr (k) near 
fcA, = 1 , and to shift 7 (k) to larger k. 

The dependence of ui on /r/T for fixed fcA, is shown in 
Fig. [21 and it is seen from examination of uJr that there 
are small- and large-/r regimes connected by a transition 
region. This transition is caused by the change in shape 
of the distribution function from Maxwellian to Fermi- 
Dirac. The scaling of the distribution function with /i is 




Figure 1. Dependence of the real (a) and imaginary (b) parts 
of the frequency on wavenumber in units of the inverse gener¬ 
alized screening length, fcA« (^), with /r/T fixed. Black (solid): 
fj,/T = —5, blue (large dashes): /r/T = 0, orange (dashes): 
/r/T = 5, red (dots): ^jl/T = 10. 


thus accounted for by the change in A*. 

Since the primary goal of this paper is to understand 
the dependence of Landau damping on the level of de¬ 
generacy, a quantity of particular interest is T (/x) = 
limfe_>oo 7 (fc, m) /wr(fc,M)- We find that the quantity T 
is approximately constant for small /x/T and then turns 
downwards to decrease inversely with degeneracy for 
larger ^/T approaching T (/x) ft! 1.3 

We have confirmed that the results obtained with our 
general numerical method agree well with the previously 
published analytical results in various limiting cases. For 
example, our calculations reproduce the analytical ap¬ 
proximations obtained by Melrose and Mushtaq for 
the dependence of both cOr and j on k and /x in the 
long-wavelength A:A* <C 1 limit for arbitrary degeneracy. 
Also, in the degenerate limit (/x/T —> 00 ) our present 
results agree well with the solutions for uJr derived for 
lon g an d moderate wavelengths (small- and moderate k) 
by ia,[ii using a completely degenerate top-hat equi¬ 
librium distribution with Fq (v > vp) = 0. For larger 
wavenumbers, Ref. 0 claimed that there is a critical 
maximum wavenumber above which no solutions exist in 
the fully degenerate case, due to the non-analiticity of 
the distribution function. Our present study (conducted 
for plasmas with arbitrarily large, but finite degeneracy) 
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Figure 2. Dependence of cdr (a) and the relative damping rate 
uJr/^ (b) on level of degeneracy with feA* fixed. Red (dotted): 
kX, = 0.5, orange (dashed): fcA* = 1 blue (large dashes): 
kXt = 3, black (solid): kXt = 8. 

cannot confirm this result for the parameter regimes we 
have explored. 

Discussion and Conclusions. In conclusion, in this 
paper we analyzed the complex frequency spectrum of 
linear longitudinal electrostatic waves in an arbitrarily 
degenerate electron plasma with a stationary neutraliz¬ 
ing positive background. Using an appropriate analyti¬ 
cal form for the dispersion relation from quantum kinetic 
theory, we developed a numerical procedure to solve for 
the real part of the frequency ujr and the Landau damp¬ 


ing rate 7 as functions of the wavenumber k and the level 
of degeneracy fi/T. We found that there are two ways 
in which degeneracy affects the dispersion relation: ( 1 ) 
through the change in the shape of the distribution func¬ 
tion from Maxwellian to a Fermi-Dirac profile, and (2) 
the increased characteristic velocity as a function of p/T. 
Above a certain level of degeneracy, the distribution func¬ 
tion assumes an approximately constant shape; beyond 
this, the real part of the frequency only depends on de¬ 
generacy through the scaling of the characteristic wave¬ 
length A*. By accounting for this scaling we isolated the 
effect of the shape of the Fermi-Dirac function on the 
complex frequency w = Wr — 17 . 

We showed that for wavelengths shorter than A, 
both LOr and 7 grow linearly with k for all levels of degen¬ 
eracy with the relative damping rate "f/tOr decreasing as 
1.3 (p/F) independent of k as p/T —>■ 00 . This means 
that electrostatic waves in such systems propagate but 
have finite damping rates. Since any real system occurs 
at finite temperature, the choice of whether to account 
for the presence of electrostatic waves, and whether to 
account for their damping, should be informed by the 
specific temperature and density of a degenerate system. 
This effect will be of importance not just for pure elec¬ 
trostatic waves, but for other plasma waves in which col¬ 
lisionless kinetic damping occurs. 

A thorough understanding of Landau damping in the 
linear regime is important before moving on to nonlinear 
effects. The realization that nonlinear particle trapping 
is suppressed in certain parameter regimes in quantum 
models gives additional impetus to study the linear 
regime of Landau damping, in contrast to classical plas¬ 
mas where the damping more quickly gives way to non¬ 
linearity. Additionally, the method used in this study can 
be easily extended to further studies of linear waves and 
instabilities in quantum Fermi-Dirac plasmas. Specifi¬ 
cally, the inclusion of additional populations of electrons 
or of mobile ions can be pursued in order to understand 
streaming instabilities and ion-acoustic waves in quan¬ 
tum plasmas. A study of these instabilities will be pre¬ 
sented in a future publication. 
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